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SYNTHESIS OF FEEDBACK SYSTEMS WITH LARGE PLANT IGNORANCE FOR 

i 

PRESCRIBED TIME DOMAIN TOLERANCES 


ABSTRACT 

There is given a minimum-phase plant transfer function, 
with prescribed bounds on its parameter values. The plant 
is imbedded in a two-degree-of freedom feedback system, 
which is to be designed such that the system time response 
to a deterministic input lies within specified boundaries. 
Subject to the above, the design should be such as to min- 
imize the effect of sensor noise at the input to the plant. 
This report presents a design procesure for this purpose, 
based on frequency response concepts. The time-domain 
tolerances are translated into equivalent frequency response 
tolerances. The latter lead to bounds on the loop transmis- 
sion function L(jaj) , in the form of continuous curves on 
the Nichols chart. The properties of L(jto) which satisfy 
these bounds with minimum effect of sensor noise, are de- 
rived. The design procedure is very transparent, providing 
the designer with the insight to make any necessary trade- 
offs, at every step in the design process. The same design 
philosophy may be used to attenuate the effect of distur- 
bances on plants with parameter ignorance. 



CHAPTER I 


STATEMENT OF THE PROBLEM 


1. 1 Introduction 

This report is devoted to the following problem. There 

is a single input-output 'plant' imbedded in a linear 

' two-degree-of-freedom feedback structure'. The term 'plant' 

denotes the constrained part of the system, whose output is 

the system output. The designation 'two-degree-of-freedom 

feedback structure' indicates a system wherein the command 

input, r(t) in Fig. 1.1, and the system output c(t) , may 

be independently measured. In such a structure 1 the system 

response to command inputs, and the system sensitivity to the 

2 

plant may be, to some extent , independently controlled. The 
structure shown in Fig. 1.1 is of course only one of many 
possible canonic 1 two-degree-of-freedom feedback structures. 
The plant parameters are not known precisely. Only the 
ranges of their values are known. For example, the plant 
transfer function may be a known function of elements of the 
set X ={x 1 , x 2 ,...x } but these elements are known only to 
lie in a given closed region in n-dimensional space. Clearly, 
the 'bounded ignorance' of the plant may appear in many dif- 
ferent forms. Strictly speaking, the design technique is ap- 
plicable only to fixed parameter plants, but it is well-known 
that for engineering purposes, it is also applicable to 
1 slowly -varying ' plant parameters. It is much less known that 
the feedback is quite effective even for rather fast-varying 

3 

plants . However, this is a topic requiring considerable 
separate treatment, so it will be assumed in this report that 
the plant is fixed but there is on the designer's part ig- 
norance, in the sense previously described, of the plant 
parameters . 



F(s) 


G(s) 


P(s) 


^ 0-« N ( no j se ) 

Fig. 1.1: A canonic two-degree-of-f reedom feedback structure. 

(Darker portion is constrained part.) 


Bounds on c(t) 
^^-Bound on |c(t)| 


Fig. 1.2: Example of tolerances on unit step response, 
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The system sensitivity to the 'plant ignorance' is to be 
characterized by the resulting ignorance in the system time 
response to a specific deterministic time input. A unit step 
is chosen for the latter in this report, but any other input 
may be chosen for that purpose. The problem is to guarantee 
that the output ignorance is contained within prescribed 
bounds, for example, those shown in Fig. 1.2, for the case 
r(t) = u(t) , the unit step function. It is, of course, up 
to someone somewhere in the hierarchy, to provide these 
bounds. This report is concerned only with the means to 
satisfy such specifications. It can be shown that if the 
plant is minimum-phase, then any such specifications, no 
matter how narrow the tolerances, may be approached as closely 
as desired. From this, it follows that it is easy to over- 
design in minimum-phase systems. The 'price' paid in general 
is in the large 'bandwidth' of the loop-transmission function, 
with transfer function L(s) / which in turn opens wider the 
'window' to the noise in the feedback return path, lumped as 
sensor noise N in Fig. 1.1 . In the high-frequency range 
where | L ( j oi ) | < 1 , but where |P(jto)l < | L ( joj) | , the 
noise is amplified by | L ( j w ) /P ( j w ) | , which tends to be very 
large over a large bandwidth, in such systems. The amplified 
noise tends to saturate the output stages of G(s) or input 
stages of P(s) . This problem has been emphasized before^ . 

1 . 2 Optimization Criteria 

This report copes with the above problem by taking the 
response bounds as inviolate, but attempting to satisfy them 
with an L(jw) whose magnitude as a function of frequency is 
decreased as fast as possible. This decreases the effect of 
the sensor noise, if it is strong in the relatively higher 
frequency range. Another approach would be to minimize an 
index into which enter both the spread in the response and 
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the effect of the noise. Statistical methods have been very 

useful for such indices, for many problems in which the plant 

parameters are precisely known. Attempts to do the same for 

the present problem have been unsuccessful because of the 

need to obtain expectations over P of expressions like 

5 

P/ (1+PM) . The practice has been to neglect the ignorance 

of P in the denominator, and replace it there by some 
nominal P 0 , but' this is obviously a very poor approxima- 
tion for plants with significant ignorance bounds. The op- 
timization criterion used in this report is discussed 
further in Section 3.1 . 


1.3 Previous Work 


There have been two principal approaches to the problem 

considered in this report, under the constraint of linear, 

fixed compensation functions G,F in Fig. 1.1 . (This 

constraint eliminates the class of nonlinear feedback, or 

nonlinear 'adaptive' systems.) One approach is that of the 

6 7 8 

dominant pole-zero method ' ' which has been reasonably well 
documented. The other approach has been by the frequency- 
response method, of which this report is the most recent con- 
tribution. The philosophy c£ the frequency-response method 
is to translate the system time-domain bounds into equivalent 
bounds on the magnitude (and/or phase) of system transfer 
function frequency response. This first step remains to this 
day nonrigorous and approximate. The second step is to 
translate these frequency response bounds into equivalent 
ones on the loop transmission function, L(s) = GP of Fig. 

1.1 . In previous work the method used was based on the 
9,10,11 

equation 


T ( j co) 

T o (j«) 


1 + L 0 ( jco) 


P 0 

+ L 0 ( j-w) 

P 


( 1 . 1 ) 


where P Q is the nominal plant function at which the system 
transfer function has its nominal value, T 0 . When the 
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plant function is P , the corresponding system function is 
T . L 0 is the nominal value of the LTF, i.e., L 0 = GP 0 . 
From Equation (1.1) , given the range of , it is 

straightforward to find the bounds on L 0 (jco) , at any co , 
such that the bounds on T Q /T are not violated. These 
bounds on L 0 (jco) are obtained at as many discrete fre- 
quencies as one may wish. The third step is to choose a 
satisfactory L 0 (jco) , such that the above bounds on L Q 
are satisfied. This determines G = L 0 /P 0 , and then F 
is chosen from the relation 

T o = F L # /(l+L 0 ) . 

This design technique has remained virtually unchanged 
since 1959 when it was first presented. It has three 
distinct shortcomings: (1) the nonrigorous, approximate 

nature of step 1, in which time-domain bounds are translated 
into equivalent frequency-response bounds. (2) the need in 
the second step to make at the outset a pairing of a nominal 
response function with a nominal plant function. An ex- 
perienced designer can usually make a reasonable good 
pairing. Normally, the slowest admissible response function 
corresponds approximately to the slowest, smallest magnitude 
plant. However, even the very experienced designer has very 
small probability of a priori choosing a pairing which 
proves to be optimal. And the inexperienced designer may 
initially make a very poor pairing and be forced to a 
second round. This specific shortcoming disappears in the 
revised method presented in this report, by eliminating en- 
tirely the need for pairing. Finally, the third shortcoming 
is in the lack of guide-lines for the third step, i.e., 
finding the optimum L Q (jco) which does satisfy the set of 
bounds derived in step 2. This shortcoming is very signi- 
ficantly reduced in the present work, in which there are 
presented many of the properties of the ideal L(jco) . 



CHAPTER 2 


DERIVATION OF BOUNDS ON L(joj) 

2. 1 Templates of P(jio) 

It is assumed that the translation of time-domain into 
frequency-domain specifications has been accomplished and 
the latter are of the form shown in Fig. 2.1. For example, 

at w =5 rps , | C ( j 5 ) | is permitted to be anywhere between 

-2.5 and -13.5 db ; | C ( j 10 ) ) between - 7.0 and - 27.0 db 

etc. Since C = RT = RFL/(1+L) and there is no ignorance 

of F, G , 

AAn|c(jw)| = M n | T ( j a; ) j = A£ n 
with 

L = GP (see Fig. 2.2) ' 

Ail nL = Ail nP . (22) 

A specific value of frequency is chosen; say 
u = Wjrps . The values of P(jo)j) over the range of plant 
parameters are calculated and the bounds obtained. The 
procedure is illustrated for the case 

k ci 

P(s) = s (s+a)' ; 1 i k < 10 ; 1 < a < 10 . (2.3) 

This is done on the plane of Jin L (joi) = £, n | L | + j Arg L 

the abscissa in degrees and the ordinate in decibels (the 
Nichols chart). Thus, at u = 2 rps, P (2 j ) lies within 
the boundaries given by ABCD in Fig. 2.3. Since 
ilnL=ilnG+£nP , the pattern outlined by ABCD may be 
translated (but not rotated) on the Nichols chart, the 
amount of translation being given by the value of G(2j) 

For example, if a trial design of L(2j) corresponds to 


L ( j lo) 


1+L ( j ai ) 


( 2 . 1 ) 
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the template of P(2j) being translated to the new position 
A'B'C'D' in Fig. 2.3, then 

| G ( 2 j ) | = | L (2 j ) | - | P (2 j ) | = (-2.0) - (-13.0) = 11.0 DB ; 

Arg G (2 j) = Arg L(2j) - Arg P (2 j ) = (-60°) - (-153.4°) = 93.4° . 

2 . 2 Bounds on L(jo) in the Nichols Chart 

The templates of P(jo) are manipulated to find the 
position of L(jw) which results in the specifications of 
Fig. 2.1 on Jln|T(jo)| being satisfied. Taking the to = 2 
template, one tries, for example, positioning it, as shown in 
Fig. 2.3, at A'B'C'D' . Contours of constant it n| L/ <1+L) j 
are available on the Nichols chart. Using these contours it 
is seen that the maximum change in £n|L/(l+L)| which, from 
Equation 2.1, is the maximum change in n | T | is, in this 
case, very closely (-0.49) - (-5.7)=5.2 DB, the maximum being 
at the point C' , the minimum at the point A' . The speci- 
fications of Fig. 2.1 tolerate a change of 6.5 DB, at to = 2 
so | L ( j 2 ) | is in this case more than satisfactory. One may 
shift the template lower on the Nichols chart, until the 
bounds on M n | T | correspond to 6.5 D.B. This is achieved 
when the lower left corner of the template is at A' ' in Fig. 

2.3 . The template corners are at A''B''C''D'' and the ex- 
treme values of &n/L/(l+L)| are at C'' (-0.7 DB) , A’’ 

(-7.2 DB) . If L(2j) for condition A is chosen to be 

“4.2 DB /- _60 0 , then it is guaranteed that A£. n | T ( j ) | < 6.5 DB, 
over the entire range of plant parameter values. If 
Arg. L^(2j) = - 60° , then -4.2 DB is the smallest magnitude of 
L a ( 2 j ) which satisfies the 3.2 DB specification for M n | T | . 
Any larger magnitude is satisfactory but represents, of course, 
overdesign at that frequency. 

The manipulation of the w = 2 template may be repeated 
along a new vertical line, and a corresponding new minimum of 
I L A ( 2 j ) | found. Sufficient points are obtained in this man- 
ner to permit drawing a continuous curve of the bound on 



L A (2j) , as shown in Fig. 2.3 . The entire process may be 
repeated at other frequencies. Fig. 2.4 shows the outline 
of the templates and the resulting bounds on L^Cjw) , for 
a number of frequencies. The template outlines are drawn 
at any available convenient areas on the Nichols chart, 
since they are, in any case, translated later. In each 
case the permissible region is to the right of the curve. 

It is important to note that although condition A (a=k=l) 
was chosen to generate the curves, any other set of values 
of a, k could have been chosen. However, once condition 
A was chosen at any one value of to , this same condition 
A must be used for all the contours. If the bounds on L A 
are satisfied, then automatically those on all other sets 
of parameter values, as encompassed by the templates, are 
also satisfied. This means that AJL n | T ( j oj ) | will not ex- 
ceed the specifications of Fig. 2.1 . 

Since P(s) is infinite at to = 0 , the system is 
Type 1, and the zero frequency specification on A£n|T| 
(assuming it is zero) can be satisfied with any finite 
value for lim sL(s) . In practice there will be a re- 

c ->- 0 

quirement on the velocity constant which will give a lower 

bound on lim sL(s). 
s-»-0 


2 . 3 System Response to Disturbances 


The system response to command inputs r (t) (of Fig. 
2.2) is not the only response function of interest. There 
are, in most systems, also disturbances to be considered. 
The disturbance response (for D in Fig. 2.2), is 


C D (s) 


D(s) 
1+L ( s ) 


(2 


It is necessary, of course, to choose L(s) so that the 
disturbances are properly attenuated, and the technique of 
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this report lends itself very readily for this purpose, as 
detailed in Section 7 ,/2 . For the present, only one aspect 

of this disturbance response problem will be considered, 
namely, that there will generally be a constraint on the 
damping factor of the pole pair nearest the axis. This 

damping factor can be related to the peaking in 

L ( jw) 

1+L ( jco) 

Thus, using the single complex pole pair as a model, a peak of 
8 db corresponds closely to a damping factor £ = 0.2 ? 

2.7 db to £=0.4, etc. The usual constraint on the damping 
factor can therefore be translated into a constraint on the 
peak value of |l/(1+l) | . Suppose this happens to be 2 db, 
in the present example. If so, the contours in Fig. 2.4 must 
be modified, as they now permit peaking greater than 2 db . 

The required modifications are shown by the dashed lines. The 
parts of the contours rendered invalid because of the above re- 
quirement are to the left of the dashed lines. Note that a 
portion of the |t| = 2 db locus is common to all the contours. 


2.4 The single high-frequency boundary 


At large frequencies, w > > a =10 

max 

\ 

P (ju) - — 

(jco) 2 

and the template becomes a vertical line of length I Aka I 

1 1 max ' 

which is 40 db. At these frequencies, from Fig. 2.1, it is 

seen than A£n|T| is allowed to be more than 40 db . Hence, 

the only constraint on L ( j co > , for the higher frequency range, 

is that |L/{1+L) | < 2 DB , due to disturbance response. This 

gives the contour marked co h , in Fig. 2,4, which will be valid 

for all frequencies higher than u>, , with co, ~ 5a . 

h h max 
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A reasonable number of discrete frequencies must be worked 
out in tie above manner to ensure that, at the frequencies 
omitted, the design will also be satisfactory. The next step 
is to choose the optimum L(joj) which satisfies the above 
bounds. Important properties of the optimum L(joi) are 
derived in the following two chapters. 



CHAPTER 3 


PROPERTIES OF OPTIMUM L(jto) - NECESSARY AND SUFFICIENT 
CONDITIONS FOR BOUNDARIES OF NEGATIVE SLOPE 

3. 1 Definition of Optimum L(jm) and its Justification 

This chapter is devoted to the derivation of the optimum 
L(ju)) function, which is defined as follows. At very high 
frequencies 

L(jw) / (3.1) 

(jw) e 

where e is the excess of poles over zeros of L(s) . The 
optimum L(jw) is the one which satisfies the specifications, 
as obtained in Chapter 2, and whos.e K, of Equation (3.1) , 
has the minimum possible value. What is the reason for this 
definition of optimum? It lies in the effect of noise in 
the feedback loop on design practicality. 

■ In Fig. 3.1, represents the power spectrum of 

sensor and amplifier noise, referred to an indicated point 
in the loop. A typical difficult adaptive problem, i.e., 
one with relatively large ignorance or narrow response bounds, 
will require relatively large |l| over a large bandwidth. 

The principal reason is the existence of an almost vertical 
boundary line, extending over a large range, like the one 
with phase of -127°, in the specific example of Fig. 2.4 . 

This line is a boundary of permissible L(ju)) over a large 
frequency range. It extends vertically for approximately 
40 db, because the range of K of Equation (3.1) is precisely 
40 db (recall K = ka in the example of Chapter 2 with 
1 < k, a < 10) . The phase lag of ( j w) , [recall A cor- 

responds to the condition a = k = 1] , must not exceed 127° 
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Fig. 3.2: Principal factors in sensor 

noise problem 
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for w < °°h 35 5° r P s * In Fi 9* 2 * 4 ' at w = 5 , |l a | can be 
approximately — 11 DB • Let us assume L A (jto) has the maximum 
allowable lag of 127° for u > 5 up to co = 50 . Then the 
average slope of | L A | over this decade is -(127/180) 40 — 

-28.2 DB per decade. Hence, at best, | L A ( j 5 0 ) | « -11 -28.2 
=-39.2 DB . If the bounds for the frequencies between 5 and 50 
rps do not interfere, it appears feasible to have | L A ( j 5 0 ) | 

~ - 40 DB . Note, however, that if the range of k happened to 
be* say, 100 instead of 10, then the range of K would be 
60 DB and the vertical line would extend down to approximately 
-60 DB, while the boundaries at to =5,10,30 would be basically 
unchanged. In such a case, | L A | , would have to decrease at 
this 'slow' rate of -28.2 DB per decade until | L A | ~ -60 DB, 
which would require 60/28.2 = 2.1 decades beyond to = 5 , i.e., 
up to to = to x = (5) 10 * . The phase lag of L A may be 
sharply increased only after to > to v , permitting only then a 
much faster reduction of | L A | . Essentially, it is a matter 
of gain margin. If the range of K is 100 , a gain margin 
> 40 DB is required, etc. 

In Fig. 3.1, the effect of the noise at X is obtainable 
from the power spectrum of X N , denoted by 


using the relation 
^*N^RMS 

The principal factors in the noise problem can be ex- 
plained with the aid of Fig. 3.2, in which a hypothetical de- 
sign example is presented, | L A ( j 5 ) | ~ - 11 DB and decreases 

at the rate of -28.2 DB per decade up to 100 RPS at which it 
is generously assumed that the slope can be instantaneously. 


N 


1+L 


= 


N 


L/P 

1+L 


1 / 2 


7 

= { [ $ doj} 

1 1 XN 


(3.2) 


(3.3) 
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changed to -100 DB per decade. |P A (jto)| is also shown, and 

so is |G | . In the frequency range where 

l L J <<1 , ^1+L^J *' 1 , and | | < 1 for w > 40 which, 

from Fig. 3.3, is the most significant range for calculating 
(x n )rms • Thus, in this important frequency range 

*XN * *N |L/P i “ " * N I Gl 2 ■ (3-4) 

Since P is beyond the designer's control, it is most im- 
portant to reduce | L | as quickly as possible, as a function 
of frequency. Hence, the importance of having tolerances as 
large as possible for Mn|T| , and for optimizing L a (ju>) , 
in the sense indicated. It is important to note that arith- 
metic scales are used in Fig. 3.3 for finding (x^) . 

Hence the saving of even a few db or of a part of an octave 
in L ( j ol) ) has a significant effect on x N . 

The definition of optimum L(ju) , [minimum K of Equation 
3.1], is reasonable when the noise is white. But if this is 
not so, for example, if there is a sharp peaking of the noise, 
say at oo = 10 , in the example of 2.4, then it would appear 
that a more precise definition of optimum directly in terms of 
$ N etc. would be necessary, leading to emphasis in the 10 rps 
range where $ N is unusually large. Note, however, that in 
any case such a redefinition would not help matters much, as 
the bound on L(jl0) in Fig. 2.4 precludes drastic reduction 
of | L ( j 10 ) | . (In fact, it is not an unreasonable conjecture 
that the definition of optimum L being used leads to the 
minimization of | L | at every frequency, but this is only a 
conjecture.) Hence, significant reduction of | L ( j 1 0 ) | re- 
quires a significantly lower boundary on the Nichols chart of 
L(jl0), which would, of course, affect the time domain response 
for some combinations of plant parameter values. This is now 
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a matter of compromise between the time domain specifications 
and the noise response, and would have to be handled at that 
level. Of course, an optimality criterion which includes 
within it both these factors would therefore be superior in 
this respect. But, as noted in Chapter 1, there has been no 
success in applying such a criterion to this problem. 


Estimate of L(jto) . An excellent idea of how fast 
| L « ( j 0 ) ) | may be decreased is obtainable from a simple physical 
inspection of the contours of Fig. 2.4 . It is of course 

desirable that L, have as large a lag angle as possible be- 

A 13 

cause the phase is, m a general sense proportional to 

jn I j! 

-gs — ' — L . Since the boundaries tend to increase with phase lag 
d2, mo * v 

(at least on the left side which is the desirable larger phase 
lag side), there will be a compromise. If the slopes of the 
boundaries are rather gentle, as they are for those in Fig. 

2.4, then the optimum L A will probably pass through points 
near the minima on the boundaries. There is also the matter 
of consistency, of the slope of |l a | with the location of 

the contours. For example, if L A ( j 2 ) is cnosen ~ - 6.5 DB 

/ 

h-127° and it has this lag angle over an octave or so, then 
the rate of decrease of |l a | near to = 2 is approximately 
(127/180) 12 = 8.5 db per octave. Hence at to = 5 , which is 
1.32 octaves away, |l a | 'will be approximately -6.5 - 1.32(8.5) 
= -17.7 db, which is well below the boundary for to = 5 . The 
optimum will therefore probably have less than the above 

phase lag at to = 2 . It is seen that the problem is not an 


6 . 5 DB 


easy one. 


3. 2 Derivation of Optimum L(jto) for negative slope boundaries - 
necessary conditions 


A boundary curve on the Nichols chart (for example, as in 
Fig. 2.4) can be written as | F ] = f { 0 ) . It is assumed in 
this chapter that d | F | /d0 is negative over the entire range. 



This is generally the situation at larger frequencies, but 
isn't so at the lower frequencies (see Fig. 2.4) . However, 

it has been noted that it is desirable to work in the larger 
phase lag region in which all the contours have the negative 
slope character. If the system is Type 1 or 2, then one may, 
in fact, be in this negative slope region for all, or almost 
all of the frequency range. This is not true for Type 0 
systems where Arg L(0) = 0°, where L(jw) must originate 
on the zero degreee vertical line and then, perforce, must 
cross a region where the contours tend to have positive 
slope. Nevertheless, we shall first assume the negative 
slope condition, because it is then possible to determine 
necessary and sufficient conditions for the ideal L (de- 
noted as Lj) and prove that Lj exists and is unique. 

The more general case of boundaries of general V shape, is 
considered in Chapter 4. 

The derivation of the optimum L(jw) for the negative 
slope case, is accomplished in the following steps. 

(1) Derivation of necessary conditions. 

(2) Uniqueness of the function which satisfies the 

necessary conditions of Step 1, on the assumption 

of fixed L(0) for Type 0 system, or fixed 

lim sL(s) for Type 1 system, etc. This is done in 
s~*"0 

Sec. 3.3 . 

(3) Proof that the optimum, if it exists, will have its 

L(0) [or lim sL(s) etc.] at its minimum per- 
s+Q 

missible value (Sec. 3.4). 

(4) Proof of existence of optimum (Sec 3.5). 

Taken together, the above provide necessary and sufficient 
conditions for the optimum L(jw) . 
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Step 1 . Derivation of necessary conditions 

Necessary conditions are derived by proving that any 
postulated Lj (jw) , which does not satisfy these conditions, 
can be improved by using instead 

L 2 ( joo) = F ( ju) Lj ( joj) 

which is "better" than Lj (joo) . For this purpose the three 
F (F,. , F_, F 0 ) building blocks in Fig. 3.4a, b, c are useful. 
It can be assumed that the phase in Fig. 3.4a and the magni- 
tude characteristics in Fig. 3.4 b,c are first chosen and then 

the missing characteristic is obtained from the Bode integral 
14 15 . 

relations ' . The important properties of these building 

blocks are indicated by the dashed vertical lines in Fig. 3.4 , 
viz., for all oo to the left of XX' in Fig. 3.4a , / F^ is 
zero, while I f jJdb 51 0 etc * In Fig. 3.4a, note 15 that 
A = (2 /tt) times the shaded area of Fig. 3.4a . 

Let M x be the locus, on the Nichols chart, of the 
maximum |l/( 1+L)| allowed due to the peaking specifications 
on the disturbance response. Let V be the vertical line 
tangent to , corresponding to the maximum phase lag al- 

lowable until the required gain margin is obtained to take 
care of | AK | (of Equation (3.1)). For example, V extends 
from 5 to - 36 DB in Fig. 2.4. Let be the smallest 

frequency for which V constitutes the right-hand part of the 
boundary of permissible L(joi^), i.e., for all co > co^ the 
boundaries have no portions to the right of V (u>^ ~ 50 in 
Fig. 2.4). It is generally the lowest frequency for which the 
permissible AHn|T| exceeds the height of the P template. 
This usually occurs at a frequency sufficiently large such 
than the P template is almost a vertical line. 
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Necessary condition 1 (for boundaries of any slope) 

For all w > , the optimum L(joo) will have no portion 

to the right of V . 


Proof. 


Suppose an L x (jco) is postulated which does not satisfy 
Condition 1, but for some frequency range < co < w 2 , 
lies to the right of V , as in Fig. 3.5 . This means that 

Let 


Lj for this range has less phase lag than that allowed. 


with 
the line YY' . 


L " F A L i 


0 ), 


on the line XX' in Fig. 3.4a, 


to , 


on 


The arrows in Fig. 3.5 indicate the effect of 


f a on 


L . For all a) < to 


at the same phase, 


so that Lj certainly lies on or above the boundaries if Lj 
does (see Fig. 3.5, where the dots denote L 3 


positions and 

the arrows denote the general position of L = F A L X ) . For 
to, < to < oj 2 , Jh -C Jl x , but this is tolerable. The level J 


i - 


of the building block of Fig. 3.4a need only be adjusted such 
that L does not cross the \ertical line V . J is finite, 
even though it may be small. At high frequencies |l| < [ L x | 

by a finite amount A = (2 /tt) times the shaded area of Fig. 
3.4a . For all to > io 2 J L = J L } , so the constraints on L 
for to > to are still satisfied. Note that this necessary 
condition applies, whatever the slope of the boundaries, as 
the proof is valid for boundaries of any slope. 


Necessary condition 2 (for boundaries of any slope) 

The optimum L(jto) must lie on the boundaries in the 
frequency range immediately below to^ . 

Proof 

Again, let L not satify this condition, in that for 
the range < w < it lies above its boundaries (see 

Fig. 3.6) . Again, let L = F,L using building block A 

of Fig. 3.4a for F^ , with on XX ' , on YY ' . 
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The effect of F. is indicated by the arrows in Fig. 3.6 . 

For all a) < (jo , L certainly lies above the boundaries, if 
L did. At high frequencies |l| < | L : | by the amount A 

= 2/tt times the shaded area of Fig. 3.4a . Only in the range 

a) < w < , does L tend to be closer to the boundaries than 

x - - h 

L , but the latter has a 'surplus' area in precisely this 
range. J may be adjusted until L lies precisely on a 
boundary at one or more < w < . Again, the proof is 

valid whatever the slopes of the boundaries. 

Necessary condition 3. (only for boundaries of negative slope) 

In the frequency range immediately above u> = 0 , the op- 
timum L must lie on the boundaries. 

Proof . 

Assume L x does not lie on the boundaries for 0 < to < w , 

(Fig. 3.7) . Let L = F b L i • F B of tile shown in Fig. 3.4b, 

with a) on XX' . The effect is shown by the arrows in Fig. 

3.7 . For all to > u 1( the new L lies above the boundaries, 
even if L x lay exactly on the boundaries for to > to 4 
Hence, the techniques used for conditions 1, 2 may now be used to 
improve matters. Essentially, F B permits the shifting of a 
"surplus" from the low frequency region, to a more convenient 
high-frequency region, where the techniques of conditions 1,2 may 
be used. 

Necessary condition 4 (only for boundaries of negative slope ) 

The optimum L cannot have any range in which it lies above 

the boundaries. 

Proof . 

It has already been shown that the optimum L must lie on 
the boundaries above and immediately preceding , and for the 

range immediately following w = 0 . Consider an L which 




27 . 


lies above the boundaries for any range oo < oo < oo n (Fig. 
3.8) . Let L = F c L x , F c of Fig. 3.4c with oo A at XX' , 
oo R at YY' . The arrows show the effect of F in Fig. 3.8. 

For all oj < oj a and oo > oo R the effect is to move the new L 
away from the boundaries, thus creating pockets of "surplus" 
which may be exploited by the techniques of Conditions 1-3 . 

In part of the range oo A < oo < R , the effect of F c is to 
push L towards the boundaries but a surplus exists in this 
range for this purpose. 

From the above four conditions, for boundaries of nega- 
tive slope, it follows that the optimum L , if it exists, 
must lie at each frequency, on its boundary. The next step 
is to prove that such an L , if it exists, is unique, if its 
value is fixed at oo = 0 . 


3. 3 Uniqueness of L with prescribed L(0) 

Theorem . 

There cannot be two or more L(joo) functions, which lie 
on the boundaries for all w >0andv#ich are equal at oo= 0. (Proven only 
for boundaries of negative slope.) 


Proof. 


and L 


Suppose that it is not so and there are two of them, L 


B 


with 


l a (0) = l b (0) 


There are three possibilities. 


as shown in Fig. 3.9 


in which case 
for all 


J B1 


all go 
case 

■ & 

W * L B3 -*- s for t]le last possible case when L R 


I l b' 1 L A 


is for the case 

for all oo . 
, in which case 0 


e B i 8 A 


00 


B2 


B - 0 A 


I L a | for 

is for the 

for all 

crosses L, 
a A 

at one or more points. (Note that at oo = 0 , in a type zero 
system, L(0) must lie on the vertical line 0=0. In a 
type 1 system there is a lower bound on lim, sL(s) , and 
sL(0) must lie on the vertical line 0 =~*-90°.) Consider 






L. = L /L . for the three L_ cases. Magnitude and phase 
1 A Bi ^ 

sketches of I L i I DB and e i are shown in Fi 9* 3 - 10 ' for 
the three . 


Apply the relation 16 (which applies to minimum-phase 


L i> 


i / 6 i iln|L i |^j = (ft n | ( j 00 1 ) 2 - (2- n | L ± (0 ) | ) 


(3.5) 


to each L. . Since L i (0) = 1, the righthand side of (3.5) 
is positive. However, from Fig. 3.10, it is seen that when- 
ever 0- > 0 , ftn|L-| 1 0 and vice versa. Hence the left- 
hand side of Equation (3.5) must be negative. Thus the sup- 
position is untenable and there can be only one such loop 
transfer function. It is very important to note that this 
proof depends on boundaries with negative slope as in Fig . 


3.9, for only then is 


when 


-B 




The above theorem does not preclude an infinitude of 
Lj ( j u) , providing each has a different value at zero 
frequency. 


3 . 4 Optimum L has minimum permissible value at u) = 0 

\ 

The next step is to prove that the best L is the one 
with minimum magnitude at w = 0 (or for Type 1 systems the 
one with minimum value of lim sL(s) , etc) . This is done by 
means of the following relation, for minimum— phase F(s) , 


c [ft n | F ( 1 a ) | -il n | F (°°) 1 ] dm 




0doo 



(3.6) 


where 





Theorem . The optimum L must have the minimum permissible 
magnitude at to = 0 . 

Proof . 

Consider L (s) with (1) the minimum permissible magni- 
m 

tude at to = 0 and (2) which crosses the boundaries at each 

value of to . Compare it to any other Lj(s) which satisfies 

(2) but is larger than L m at to = 0 . L n ( j<o) can cross 

L (jto) either an even or an odd number of times in the range 

0 < to < to, . If they do not cross at all in this range, then 
h ^ 

the situation can be as shown in Fig. 3.11, for F = L l /L m . 
Case (a) or any even number of crossings in the range to > to h 
leads to 

! p L = IW 00 * ! - 1 ' (3 

so there is no advantage in L ; over L 2 . Case (b),, ap- 
plied to Equation (3.6) with u> c = u> h , has zero for the 
right-hand side. On the left side &n|F(°°)|< 0 . Hence 
the integrand is always positive. This is a contradiction 
in that the left side is finite and positive, the right side 
zero. Hence case (b) is impossible. We next turn to the case 
of even number of crossings in the range 0 < u> < u)^ . 

Case 1 Even number of crossings in the range 0 < to < to^ 

In this case F(jw) = Lj ( jto) /L ( jto) will appear 
qualitatively as in Fig. 3.12. Sub-case (i) or of any even 
number of crossings in the range to > to^ need not be con- 
sidered, for then (3.8) applies. Only sub-case (ii) , i.e., 
an odd number of crossings in the range u) > to^ need be 
considered. I 


l 






34 . 


Let u)^ = to of Equation (3.6) and write &n|F(°°)| = -M , 

M > 0 because, for subcase (ii) , |f(°°)| < 1 . Equation (3.6) 

can be written as follows 


(otj-otj) + (a 3 -a 4 ) + • * * + (a 2n-l -a 2n ) + V i + U 2 = 0 


(3.9) 


where 


to 


a 


2i+l 


= / 


2i + l 


(0 


2i 


in F din 


1 - — 


CO. 


h 


(3.10a) 


a 2i " 1/ 


ID. • 
2i 


in F da) 


OJ 


2i 




10 


10 , 


(3.10b) 


to. 


= / 


i n I F dt) 


W 2n Jl - JL 


to. 


(0, 


= M / 


dto 7T 


/l - 


= ^-Mto, > 0 
2 h 


to; 


(3.11a) 


(3.11b) 


Next, let lo „ = to of Equation (3.6), which becomes 

^ n c 


+ (3 3 -3 4 ) +•••+ ce 2n _ 1 -e 2n ) 


+ 


y 


2 


£ > 0 (3.12) 
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because on the right-hand side of (3.6), 6 < 0 in the range 

w = to_ < to < to, . In (3.12), 
c 2n h 


to 


e 2i+l 


« / 


2i+l 


(0 


2i 


in F dm 


!i - 


(3.13a) 


<0 


tO 


2n 


’ 2i 


(0 


= 1 / 


2i 


(0 


i n ! F d-i 


3.13b) 


2i-l 1 - 


to 


0) 


2n 


The last two terms on the left-hand side of (3.12) constitute 


0) 


- / 


2n 


% n F (°°) I dto 


1 - 


to 


TT TT 

2 2n 2 2h 


(3.13c) 


V 


co 


2n 


which can be written as ^ Ma) 2h ” G = ^ 2 “ E with e > 0 . 

\ 

Compare any pair ^ a 2i-l _a 2i+2^ Equation (3.9) 
with (^2i+l~^2i+2^ of Equation (3.12) . 


a 2i+l a 2i+2 


<0 


- 1 / 


to 


2i+l 


/ 


ID Vxl 1 - 
2i+l u 


to_ . 

2i 


1 - 


to 


2i+2 


to 


2n 


1 - 


to 


h 


to ^ 0.5 

2n j & n | F | dto 


1 - Si- 


■/I 


to 


0.5 


V 


to 


2n 


l n F dto 


1 - 


to 


to 


2n 


(3.14a) 
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The function p(m) = 


1 

[- 

1 



2 


_^2n_] 



0.5 


(3.14b) 


is monotonically decreasing for 0 < w < w 2n . Hence, 
P( w 2 , +l) < p for the °P en interval of integration of the 

first term on the right side of (3.14a) , while P ^ u 2i+l^ 

for the open interval of integration of the second term. But 
note that if p (w) is removed from these integrals, what re- 
mains are precisely the integrands of ^2i+l ' 0 2i+2 ' re ~ 
spectively, and over the same limits of integration. Hence 


a 2i+l > p (a) 2i+l ) 3 2i+l ; a 2i+2 < p {w 2i+l } 3 2i+2 

a 2i+l “ a 2i+2 > p (u) 2i+l ) (e 2i+l" 6 2i + 2 ) 


(3.15) 


(3.16) 


Note also, because of the monotonically decreasing nature of 
p(w) , that 

! > P[ > p 2 ”* >,P 2n _! > 0 • 

Combining (3.16) and (3.17), gives 

(a r a 2 ) + (a 3 - a J + -*- (a 2 n-r a 2 n ) 

> p(to i ) (6 r B 2 ) + P(w 3 ) (e 3 -^ ) + '-' p(w 2n-l ) (3 2n-r 6 2n } 

> P(^n-l )[(3 1 “ 0 2 ) + + * * * (3 2n-l _3 2n ) - 


> P (w 2 n -i ) (-P 2 +e) ' 


(3.18) 


from (3.12) . 
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But, from (3.9), the left side of (3.18), i.e., 
(a j -a 2 ) +• • • = - (yj+y 2 ) . 

Hence, 


- ( y x +y 2 ) > p (w 2n _ 1 ) (-y 2 + e) 


(3.19) 


or 


- y j-EP ( a) 2n-l^ > ^ -p ( a) 2n-l^ > 0 * (3.20) 

because 0 < p(co 2n _ 1 ) -C 1 • 

Since y 1 , e , p are all positive (3.20) is impossible, so 
the hypothesis is untenable. 

Case 2 . Odd number of crossings 

In this case F(jw) = Lj ( joo) /L 2 ( j w) will appear quali- 
tatively as in Fig. 3.13. Sub-case (i) in Fig. 3.13, or of 
any even number of crossings in the range co > oo^ , need not 
be considered, for then (3.8) applies. Only sub-case (ii) , 
i.e., an odd number of crossings in the range oo > co^ , needs 
to be shown impossible. 

Let u> c of Equation ' (3.6) be set equal to C0 2n+1 * 

Then Equation (3.6) becomes 


( a r a 2 ) + (a 3"°\> + "* +(a 2n-r a 2n ) +a 2n+l 


, TTM 

+ ~T 0J 2n+l 


= - / 


0dco 


= - a < 0 (3.21) 


w 


2n+l 


co 


- 1 


1 


00 


2n+l 
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where the or s are defined as before, in Equations (3.10a,b), 
except for the^ substitution of L0 2n+1 ^ or w h * 

Let w = <o 2n in Equation (3.6), and it becomes 


to,. 


(6 1 -6 2 ) + (6 3 -0 4 ) + ---+(6 2n _ 1 -6 2n )+^ 03 2n = / 


0dto 


= -a+B, 


(0 


2n 


to 


-1 


to 


2n 


with B > 0 , since 


to,. 


to. 


a = / 


0dio 


0dto 


to 


2n+l 


to 


, to„ 
-1 2n 


to 


-1 


to 


2n+l 


to 


2n 


because u 2n+l > w 2n and because 9 < 0 in the range 

w 2n < “ < u 2n+l * 

As in the previous development, 


(a 2i+l _a 2i+2 ) > X (aJ 2i+l } (3 2i + l -e 2i+2 } 


with 


i - 4 


X (to) = [■ 


to 


2n 


1 - 


to 


w 2n+l 


0.5 


monotonically decreasing from unity over the range 
0 < to < to 2n , i . e . , 


1 > X(to i ) > A ( to 2 ) > X ( Q) 2n-1^ > 0 * 


(3.22) 


3.23a) 


3.23b) 
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Therefore 


(a i -a 2 ) + *--+(a 2n _ 1 -a 2n )>p(a J2n _ 1 ) [(6 1 -B 2 ) + *--+(6 2n _ 1 -6 2n )] . (3 
From Equations (3.21, 3.22), this means 


ttMw 


- a - 


2n+l 


ttMu 


- a 2n+1 > pfVl 1 1- “ + 6 


2n. 


(3 


which can be manipulated into 

7T Moj r\ jr fcfl 

<“+— [l-P(“ 2n - 1 )l + T-<“ 2n+ r“ 2 n ) <- a 2n+r Bp( “2n-l) ' (: 

Each factor of each term on the left is positive, while each 
term on the right is negative. Therefore the assumption 
(ii) of Fig. (3.12) is impossible. 

So far, it has been found that the optimum L , if it 

exists, must lie at each frequency on a boundary curve and 

have the minimum permissible value at u> = 0 . It has also 

been proven that such an L , if it exists, is also unique. 

There therefore only remains to prove that such an L(s) 

, \ 
does indeed exist. 

Existence Theorem 

There exists an optimum L function in the sense 
previously described. 

Proof 

We use the analogue of the Weierstrass Theorem which 
states 18 that if J(f) is a continuous functional defined 
on the normal and compact family {f(z)} , then |j(f)| = MIN 


.24) 


.25) 


.26) 



has a solution within this family. In this problem 

J(f) = J C L ( j oi ) ] = lim| (jto) x L(j<o) | , where x is chosen so 

0J->oo 

that J is nonzero and finite; for example, if the vertical 

boundary for to > to^ is at - 8 , then x = Z0 /tt . Such a 

functional is clearly continuous for a large variety of 
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norms. There is a theorem which states that a class of 
analytic functions regular and uniformly bounded in a domain 
D is a normal family which is also compact. Uniformly 
bounded functions f(s) in D are those for which there 
exists a positive constant M , such that |L(s) | < M , 
for s in D . In this case D is the right half-plane 
and L(s) is regular in D . For any given problem it is 
always possible to select an M such that the family of 
candidate L functions is restricted to |L(ito)| < M for 
all to . [For a Type 1 system, one would deal with sL(s), 
etc.] By the principle of the maximum modulus, |L(s) | < M 
in D . Hence, all the conditions of the theorem are satis- 
fied and our functional must have a minimum. 



CHAPTER IV 


PROPERTIES OF OPTIMUM L(jco) - NECESSARY CONDITIONS FDR BOUNDARIES 
OF BOTH POSITIVE AND NEGATIVE SLOPE 


4 . 1 Introduction 

This chapter is devoted to the derivation of some 
properties of the optimum L(ju) function, for the case when 
the boundary curves can have portions with positive slope, as 
for example in Fig. 4.1. In contrast with Chapter 3, a 
complete theory has not been developed. An important necessary 
condition has been found, and exactly the same existence theorem 
as in Chapter 3 applies to this case. However, sufficiency, 
proven implicitly in Chapter 3 by proving uniqueness, has not 
been proven here for the function which satisfies the necessary 
condition. 

4 . 2 Some necessary conditions for the optimum L(s) 

The first two conditions were proven in Chapter 3, as 
necessary conditions 1,2 . 

(1) For all oj > , the optimum L(jw) must lie on the ver- 

tical line V (Fig. 4.1). 

(2) The optimum L must lie on the boundaries in the fre- 
quency range immediately below . 

(3) The optimum L must be on the boundaries in the region 
where the boundaries have negative slope. 

The third necessary condition is restricted to the case where 
L(jw) passes initially through a region of positive slope 
boundaries, followed by a region of negative slope boundaries. 

Proof 

The building block F : of Fig. 4.2 is used in the proof. 

Suppose the postulated L lies above some positive slope 
boundaries, as shown in Fig. 4.1 . Then L 2 = F L will have 
the values in the direction of the arrows, and will have a 
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for necessary condition 3 
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finite ' surplus 1 for all tu > co^ . The technique of 
necessary conditions 1,2 can be used to exploit this surplus. 
Note that for the proof to be valid it is necessary that 
co ’ = (Oj , and that the range of the building block 

of Fig. 4.2 be chosen so that it lies inside the range 
(co„, to,) of Fig. 4.1 . This is evidently always possible to 
do with the many parameters available, i.e., the various 
levels and slopes of Fj . 

Generalization . More complex building blocks are required, 
in order to generalize this result to the case when L ( j co ) 
passes through several alternating regions of positive and 
negative slopes. For example, that in Fig. 4.3 may be used 
for the following order of slopes (starting from co = 0) : 
positive, negative, positive, negative, and with a surplus 
pocket in the last negative region, as shown in Fig. 4.4 . 
Obvious extensions of Fig. 4.3 would apply for larger num- 
ber of pairs of alternating slopes with the last being nega- 
tive, and wherein a surplus pocket exists in the last 
region. 

4. 3 Necessary condition 4 

This condition applies to the case where L(jco) passes 
initially through a region of positive slope, followed by a 
region of negative slope boundaries. 

(4) The optimum L must lie on the boundaries in the region 
of positive slope. 

Building Block 

A somewhat more elaborate building block is needed for 
the proof of (4) . It is constructed from Fig. 4.5, where 

1 1, | is specified to be zero db from co = 0 to co = co^ , 
while 0 = /l is specified with the indicated properties for 
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to < co ; i . e . , 0 < 0 for <o A < to < <o x , 0 > 0 for 

03 < 03 < to„ , 0 = 0 for 03 > o 3 _ . Also 0 is continuous 

X B B 

so 0=0 at 03 = o) , 03 , 03 . The missing portions 

X xJ 

( I L I for o) > o), , and 0 for 03 < 03 ) , may be obtained 

^ 20 /I 

from the appropriate Bode integrals ' , which here simplify 

into 


0 ) 


B 


0 d 5 


03 a (0) 2 -S 2 ) U 2 -03 2 ) 0 * 5 


— 7T0 (03) 
203 (o> a -03 2 ) 


for 03 < 03 


A 


(4.1a) 


M ( 03 ) 

O o 0 5 

203 (03 2 -01 2 ) * 


for 03 > 03 a 


(4.1b) 


where 

M ( 03 ) = Z n| L ( j 03 ) I (4.2) 

It will be shown that it is possible to choose 0 of the 

form indicated in Fig. 4.5, in the range 03 A < 03 < 03 g , such 

that the complete L is as shown in Fig. 4.6 . The important 

features are: the range of positive M for to > to > to , the 

y x 

ability for log |f| to be zero at infinity, or alternatively 

to cross at 03 > <o_ from positive to negative values; and to 
z a n 

remain negative from <0 to infinity. 

The proof of the above assertions will be given by 
selecting a specific form for 0 (to) , with sufficient parameters 
to obtain these properties at prespecified frequencies. 

Selection of 0 (to) in range to^ i w i W B 

Let 0(5) = /? 2 -to^ (A + B5 + C5 2 ), for u> A < to < a3 fi , (4.3) 

with 0 (co A ) = ® (“ x ) = 6 ( U B ) = 0 • 


(4.4) 
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Hence 


0 = A + Bco + Cto 2 ; 0 = A + B<o n + Cto 2 
xx an 


(4 . 5a, b) 


Solution of (4.5a,b) gives 


B = - A 


'Vb 1 

“x“b 


c = 


A 


to to_ 
x B 


(4 . 6a, b) 


so that only A and to remain to be chosen. Consider the 
sign of M(°°) which we would like to be negative or zero. 
From (4.1b, 4.3, 4.6a,b), this requirement is 


( 0 . 


B 


A / [1 - 


( W g 

W x W B 


(0 to_, 
x B 


]d£ < 0 


(4.7) 


It is easy to see that A < 0 in order that 


0 (to “) > 0 , because 
3 


6( V ) = 


A[U x W B _(U3 x +aJ B )a) B -+(a) B ) ] A(aj x ~ ai B~ ) (ai B~^B } 


(0 to_ 

x B 


to <o„ 
x B 


With A < 0 , (4.7) gives 


to' (3to'-to') - (to ' + 3oo ' ) + 2 > 0 (4.8) 

B X o o X 

\ 

where the primes indicate normalization with respect to W A ' 
i.e., <o£ = W B / W A etc * 

If 3(0 1 = to ’ , then (4.8) requires 2 > to ' + 3to' , which 

X B jo x 

is impossible. A necessary condition for (4.8) is therefore 
*» 

3 ( 0 ' 

— ~ > 1 (4.9a) 

(ol 


and condition (4.8) can be written as 
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co to' co'+2 

x = P = -4 > B 


co. 


B 


“b - 3 “b 


For example, if co^ = 2, then p = 2/3 gives M(°°) = 0, 
while p > 2/3 gives M(°°) < 0 . Also obviously, p can 
be chosen to precisely position co (of Fig. 4.6) . 

For our purposes it is important that M(co “) > 0 . 

13 

Evaluating (4.1b) with 6 (£) of Equation (4.3), gives 
M = ip 


(0 


1 - 


A 


for co > co. 


(< 


CO- 


CO 


0 = - Ip 


(0 


- 1 for co < co. 


U 


ip = A 


co 


_ LB in 

7T 1 


co A -co 


+ B 2 2, n 


V“ 


co,+co 

A 


CO -CO 

- 2 (JL_A)] 

W x W B 


u 


co co^B, = co +co 
X B l X B 


CO CO 

, , X B. 

(co + ) 

CO ' 


(< 


CO CO 

x B 

CO CO B = CO +co + (co + ) 

T 2 V R 


X B 


CO 


u 


From the above it is found that M(co fi ) > 0 requires 


CO 

p = -* < 

W B 


1 

■ “B 


2co ' 

l n- 

B 


B 


Combining (4.9b), (4.11) gives 


co ' 2 + co ' - 2 
B B 


^4 


3w B (w B -1) 2tJ B 

l n -t— — r 

1+l “B 


(4.9b) 


.10a) 

. 10b) 
.10c) 

. lOd) 
.10e) 

4.11) 


1 


(4.12a) 
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which reduces to 


Z n 


2 1.5 (1-g) 

a+1 (2+a) 


(4.12b) 


where 



It is easily found 


0 ), 


co 


B 


that (4.12b) 


applies for 


0 < a < 1 


(4.12c) 


If Equations (4.10) are expanded at u> = co A + 6 , the 
result is M ■* KA /S l n 6 , with K > 0 . Since A < 0 , 

M > 0 at co = co + 6 . Similarly, it is found that 0 < 0 
at to = to, — 6 . At (0 = 6, 0 is negative, if 


( 1 + £> 


3-1 

l n3 


(4.13) 


where p = < 1 , 3 = — , which is always satisfied, as 

g > 1 . it is known that £oth M and 0 of Fig. 4.6 are 
continuous and the shape of the characteristics for frequencies 
in between the above calculated points, as shown in Fig. 4.6, 
has been checked by numerous computer runs. 


Application of Building Block to prove Necessary Condition 4 

\ 

In Fig. 4.7 there is sketched a set of boundaries and an 
L function with the properties under consideration. Let 
L 2 = FL with F like that in Fig. 4.6. Choose oo A > (^ , 
co < co' , and co, < oo . The values of L ? are in the direc- 
tion of the arrows in Fig. 4.7 (for the case these relations 
are satisfied with the equality sign) . Hence, clearly L 2 
is better than L . 




Generalization 


Consider the sequence of frequency ranges with a trial 
L crossing in order, boundaries of positive, negative, 
positive, negative slopes, with a surplus in the last posi- 
tive range. Part of the building block of Fig. 4.2, pro- 
perly combined with that of Fig. 4.6, could presumably be 
used for a building block. But this has not been explored. 

Summary 

It has been shown that a necessary condition for an 
optimum L which passes through a range of positive bound- 
ary slopes, followed by one of negative slopes, is that L 
must lie on the boundaries at every frequency. The exist- 
ence theorem of Chapter 3 is also valid and indicates that 
an optimum L does exist, but, of course, it need not 
cross the boundaries in ihe above manner. It may be such 
as to cross only negative slope boundaries or to cross 
alternately more than two different sign slopes. Intuition 
and experience suggest that for all cases the optimum must 
lie on the boundaries. 



CHAPTER 5 


SHAPING OF THE LOOP TRANSMISSION FUNCTION 


5.1 Introduction 


It is assumed that the boundaries of L for a suitable 
number of discrete frequencies have been obtained, as in 
Chapter 2. From the work in Chapter 3, 4 it is concluded that 
the optimum L is one which lies on the respective boundaries 
at each frequency. The objective is to achieve such an L(s) 
function. Some guide lines are presented here for this pur- 
pose. 

5. 2 Detailed Shaping of L(joo) 

The boundaries shown in Fig. 5.1 are used to illustrate 
the procedure with oo h = 8 rps , i.e., for all w > 8, L must 
lie to the right of the boundary marked V in Fig. 5.1 . It 
should be recognized that the phase at zero frequency is fixed 
at 0° , or -90°, etc., according as to whether the system is 
Type Zero, or One, etc. This may not correspond to the op- 
timum L(joo) , if there was no such constraint. However, the 
'loss' due to this constraint can be made very small, because 
it is possible, at as low a frequency as desired, to have 
L(joj) at the more favorable point, wherever it may happen to 
be. The practical limitation is the inconvenience of having 
poles and zeros very close to the origin. In Fig. 5.1, what- 
ever may be the value of L at zero rps, it is easy to ar- 
range that L(j 0.5) have any desired value. 

Consider the boundaries of Fig. 5.1 . Arbitrarily, let 
us try to achieve, say, -26 db at oo = 8 . Why is this im- 
practical? Note that at oo = 4 , approximately - 10 db is 
needed for j L | . In order that |l| decrease to - 26 db 
at oo = 8 , L must have an average slope of - 26 + 10 

= - 16 db per octave over this range of frequencies 
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(between co = 4 and to = 8) . But the average phase of L 
over this range is at best approximately (-180 - 130)/!= -155° . 

With such a phase there is associated an average slope of 

- 155/180 * 12 * - 10 db per octave, which is significantly 
less than the assumed - 16 db per octave. A value of -18 db 

for | L | at co = 8 seems more reasonable. Let us check this 

value in -the same manner. The average phase is; at best, now 

-130° (see Fig. 4.1), so the average slope of |l| is 
-130/180 x. 12 = db/octave. If | L ( j 4 ) | = -10 db, then 

using the above slope, | L ( j 8 ) | = -10 -8-j = -18-j , which is 
reasonably compatible with the assumption of “18 db at co = 8 . 


It is seen that the value of -10 db at co = 4 is a good 
starting point. It would be possible to use | L ( j 4 ) | = -11 db , 


providing / L ( j 4 ) = -110 


But then the average phase from 


a) = 4 to co = 8 would only be 


-110-130 


= -120° , with an as- 


-120 

sociated average slope of — ^ g x 12 = -8db /octave, so that 
| L ( j 8 ) | ~ -11 -8 = -19 , which is so close to -18-| , as to make 
little difference. What of the other boundaries? What points 
on co = 2 , 1 etc. diould we aim for? On the co = 2 boundary, 
the -130° line is clearly the best, because the boundary is 
almost flat there. In any case, the average slope for |l| 
near co = 2 is certainly such that if | L ( j 4 ) | = -10 db , then 
| L ( j 2 ) | will be greater than -8 db, i.e., L(j2) might as well 

be on V . Similarly, for L ( j ) As for L(j 0.5) , we need 

only make certain that b(j 0.5) is close to -110° or so. 

Working backwards from | L ( j 4 ) | = - 10 db, this means that 

[ L ( j 2 ) j will be approximately -1 db , [ L ( j 1 ) | approximately 

7 db , and | L ( j 0.5) | approximately 15 db , providing we ar- 
range to have approximately -130° phase for L from co = 8 to 
1 , and about -110° phase at co = 0 . 5 . 


How can the above phase values be simply obtained? An 
average phase of -130° may be obtained by alternating a lag 
corner frequency (lacf ) with a lead corner frequency (lecf) . 
Let 1+a be the number of octaves under consideration, with a 
slope of -6 db/octave over one octave and a slope of 
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-12 db /octave over a octaves. Then the average phase lag is 

C/l) * - (6 +^ . 12) 180 , (5.1) 

L_ av. 1+a 

which is to be -130° in the above example. Solving (5.1) 
for a gives a = 0.8. Thus, if we allow a slope of 
-12 db/octave for one octave duration, we should allow a 
slope of -6 db/octave for 1/0.8 = 1.25 octaves. 

The above results are applied to the example of Fig. 5.1, 

as follows. About -110° is desired at a) = 0 . 5 . This is a 

Type 1 system, so assign a lacf at , a lecf at 2u> x , 

so that the asymptotic slope of -12 db/octave is over one 

1 25 

octave, a lacf at 2 ' (2u)j) , so that the asymptotic slope of 

- 6 db/octabe is over 1.25 octaves, etc. Try various values of 
a) , until is found such that the net result of an in- 

finite series of the above leads to -110° phase at to = 0 . 5 . 

1 25 

To simplify the numbers, we used 2.5 instead of 2" = 2.37 

and set w = 1 , which is a conservative choice. In ac- 
cordance with the above procedure, the lacf at w = 1 is fol- 
lowed by a lecf at oj = 2 , a lacf at u) = (2.5)2 = 5 , and a 
lecf at (2) (5) = 10 rps . This procedure is halted at 10 rps , 
because near w = 10 we are permitted (see Fig. 5.1) a gradual 
decrease in phase. The leyel of |l| is moved vertically un- 
til | L ( j 4 ) | = -10 db. The phase lag of the design, so far, is 
Curve A in Fig. 5.2 . 

The next lacf will be set at a somewhat higher frequency 
than would result from the above formula, because we wish to 
follow it with another two lacfs, in order to have a final 
asymptote of -24 db/octaves, corresponding to an excess of 4 
poles over zeros for L(s) . The situation at w = 10 is 
examined. At present /L ( jlO) = -114°, and from Fig. 5.1, it 
can be ~ - 140°, if | L ( j 10 ) | £ - 26 db . Hence 26° more lag 
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is tolerable at to = 10 . A lacf at co = 30 contributes 
- 18.5° at a) = 10 , leaving 7.5° . The latter permits a 
complex pole pair at 4 rps, if its damping factor is 0.2 . 

The phase lag due to all poles and zeros, excepting the last 
complex pole pair, is given by curve B in Fig. 5.2, while 
the total phase lag is given by curve C . The resulting L 
is also sketched in Fig. 5.1 . The boundaries are slightly 
violated for 6 < oo < j 10 . More significantly, for oo > 12 , 

L could have more phase lag , which would permit its faster 
reduction. 

22 

The following relation gives one an idea of the magni- 
tude reduction available by increasing the phase lag: 

} 

oo j 

- / 0d £n w )= - [Hn|L(0)| - n | L («*») | ] . (5.2) 

* o t 

i 

Compare two designs, in which |l( 0) | is the same in both, 
but in one region 0 2 < 0 X by, say, an average of 20 ° over one 
octave. The difference between the two, for the left side of 
(5.2) is §( 5 ^ 3 -) fcn 2 = .154 nepers = 1.34 db . Equation (5.2) 
is useful for estimating whether, in any particular problem, 
it is justifiable to seek to improve a tentative design. Thus, 
in Fig. 5.1, at 00 = 0.5, 1, 2 the phase lag could be in- 
creased by 6 °, 13°, 6 °, respectively, so that the left side of 

\ 

(5.2) is increased in magnitude by about 1 . 7 db . However, 
most of the resulting increase in the difference between 
l n| L (0) | and £ n| L (°°) | would be due to an increase in | L (0 ) | . 
rather than a decrease in |l(°°) | . The reason is that j L ( j 4 ) | 

cannot be decreased, and since the increase in phase lag is 
in the frequency range less than 4 rps, the decrease in slope 
of £n|L(joo) | will occur primarily there. For this same 
reason, the increase in phase lag possible in the range 
00 > 12 approximately, is more appealing, for it occurs in the 



region to > 4 , so most of the effect will be in the higher 
frequency range. 
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5 . 3 Loop Shaping by Working Backwards from L(jio^ ) 


A great deal of additional insight into the loop-shaping 
problem may be obtained by considering the following example 
in detail. Fig. 5.3 presents the contours and for the present 
purpose it doesn't matter how they were obtained. The boundary 
is a vertical line for all to > 100 rps . We are interested 

, ' \ ( ,,-v 

that J’D'Cjcii ) | be as small as possible. Let us try, at random, 

|L(jl00) | = -32 db . Consider the to = 10 boundary, and try 

L(jlQ) = -10 db /- 110 ° (point Q in Fig. 5.3). Then, from 10 . 

i i n i n+i 4 n 4 n 'XCtT' 

to 100 rps, the average slope of | L | ~ - ^ — = -.-2- . 6 -7- 

db/decade, which is more than the -32+10 = 22 db needed. Hence, 
considering for the momeny only the to = 10 boundary, -32 db for 
| L ( j 10 0 ) | is achievable. In fact, | L ( j 10 0 ) | = -36.7 db is 

achievable on this basis; perhaps an even smaller value by ter- 
minating on a lower portion of the to = 10 boundary. However, 
let us be content with the point Qj on the to = 10 boundary 
(because, in any case, as will soon be seen, other boundaries 
will not permit even this value to be attainable) . 


Consider the to = 7 boundary in Fig. 5.3 and see whether 
it is possible to reach this boundary from the desirable point 
Q x , (Q being assumed the best from the viewpoint of attaining 
| L ( j 10 0 ) | min.) . The range from to = 7 to to = 10 covers 
0.5 2 octaves. Try point (0 y'-llO 0 ) on the to = 7 boundary. 

It is impossible to achieve this, because the average slope 
magnitude is only 


110 

180 


12 = 7-j db per octave, operating over 


0.52 octaves, which gives much less than the required 10 db 
from Q to Rj . In this way, one finds it necessary to move 
down to point R, (-6.6 db /-85°) , for then the average slope 

i r- *i «’> ^ 1 " 


of - [ 


110+85 


(■ 


12 


) 3 : 


2 — '180" 7 J 6*5 db /octave , over 0.52 octaves, gives 
the magnitude change of 3.4 db required between Qj and R 2 
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Next, consider the 4 rps boundary, and the transition 
from 4 to 7 rps, a range of 0.81 octaves. It is impossible 
to reach point R 2 from any point on the 4 rps boundary, as 
the smallest magnitude difference is more than 8 db, and the 
maximum average slope is less than 8 db/octave. It is there- 
fore necessary to move up higher on the 7 rps boundary, which 
must therefore make it impossible to operate at Qj , at 
10 rps. The best one can do is to find the point closest to 
R 2 which permits transition to the 4 rps boundary, because 
this will then give the point closest to the unattainable 

0 . It is found that one must go all the way up to R 

i 

(* " 2.3 /-10 2°) , with the corresponding point at S on the 

4 rps boundary, and the point, Q on the 10 rps boundary, 

- 35.5 (point W) on the -140° line, which is the 100 rps 
boundary . 

Next, consider the transition from the 2 rps boundary 
to point S . The point U (12 /-110° ) is suitable, as the 
average slope of - (— 1 - ^ — = - 8 db per octave, over one 
octave gives precisely the 8 db difference between U and S . 

In this manner a compatible gain-phase trajectory for 
L(jio) has been obtained, and which is a good approximation 
to the optimum. It is easier now to work by hand calculation, 
as it is known what to aim at. Also, it provides a good ini- 
tial input trajectory for a computer program for computing L 
and finding its pole-zero pattern. 

5 . 4 Design of Compensation, in the Loop 

It is hardly necessary to point out that the loop trans- 
mission, L , can always be written in the form L = MP , with 
M the compensation which must be inserted by the designer 
(in Fig. 1.1 , M = g) . Some set of plant parameters was chosen 
to obtain the boundaries of L . This resulting plant trans- 
fer function, at that set of plant parameters, could be denoted 



by P . If so, the resulting L(s) should be denoted by L 
n 

with M (s ) = L (s)/P (s) . M (s ) will be independent of the 

n n 

particular choice of nominal plant parameters, i.e., of P n . 
Since both L n (s) and P n (s) represent physically realizable 
transfer functions, so does M{s) . 



CHAPTER 6 


COMPLETION OF THE DESIGN 


6 . 1 Design of Prefilter F 

The work in the previous chapters, concentrating on the 
design of L(s) , guarantees that despite the specified 
bounds of plant ignorance, the change in T(jo) remains with- 
in the allowable range. Thus, if T(joi ) is allowed to 
change from -3 db + 2 db , i.e., by a total of 5 db , then 
we can be certain only that the change in | T ( j to ) | is by 
5 db . However, the range of 

L ( jtOj ) 

1+L ( j to A ) 

may be from + 2 to + 7 db , or from any value K to (K + 5 ) db . 
It is therefore necessary to make use of the prefilter F of 
Fig. 1.1 . Specifically, | F ( jto ) j must be set at - 5 db , such 
that |FL/(1 +L) | w ranges from - 3 to + 2 db . 

The above, and the problems which may arise, are illus- 
trated by the following numerical example. 

6 . 2 Design Example 

\ 

The plant transfer function is given by 

k 

P(s) = (6.1) 

s ( s 2 +S s + P ) 

P P 

with 1 < k < 1000 , and the complex pole pair of P ranging 
anywhere in the rectangle M, N, Q, U of Fig. 6.1a . The 
performance specifications require that the step response, 
characterizable by a dominant complex pole pair, be such that 



preceding page blank 


not filmed 



S plane S plane 

\ 

ig. 6.1a: Range of complex pole pair Fig. 6.1b: Acceptable range of 

Q f p ( s ) complex pole pair of 

second-order model 
of T (s ) 
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this complex pole pair lies within the region ABDEFG in the 
complex plane of Fig. 6.1b . This problem has been treated 

g 

in the literature by the Dominant Poles Method. 

The first step was to translate the dominant-poles spe- 
cifications cf Fig. 6.1 into equivalent ones on | T ( jco ) | , 
the magnitude of the system transfer function frequency re- 
sponse. This step is at present handled, in general, as 
follows. Any simple model, such as a second-order one is 
used, and the range of model parameters found which satisfies 
the domain specifications. In this example this is already 
available, in Fig. 6.1b . The next step is to find the re- 
sulting range of variation of J T ( j co ) | . These are shown in 
Fig. 6.2a with the labels corresponding to those in Fig. 

6.1b . [The reasoning is that magnitude control implicitly 
includes phase control; at least, it is worth trying. If 
not, then phase specifications can also be included and the 
design technique remains basically the same — only the 
boundaries of permissible L will be affected.] Returning 
to Fig. 6.2a, it is important to broaden the permissible 
range of | T ( jco ) j as much as possible, for clearly, this 
permits design by an L(j<o) of smaller bandwidth. To 
broaden the permissible range of | T ( j co ) | , one may proceed 

by trial and error, trying | T ( j co ) | functions which pro- 
gressively decrease faster as functions of frequency, until 
the time response specifications are intolerably violated. 
Similarly, one tries j T ( jco ) | functions which progressively 
decrease more slowly versus to , until again the time-domain 
specifications are intolerably violated. In this way, the 
significantly larger bounds shown in Fig. 6.2b were obtained. 

The procedure of the previous chapters was then fol- 
lowed. Some of the plant templates are shown in Fig. 6.3 . 
Note that for 2 < to < 10 , they extend to infinity, because 
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Fig. 6.2a: 


Frequency responses of |T(jw)| corresponding to 
Fig . 6.1b 
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Fig. 6.2b: Larger |T(jw)| bounds obtained by cut and try 























this range is included in the plant pole range variation 
(see Fig. 6.1a). The resulting boundaries on L and 
several alternative designs of L are shown in Fig. 6.4 . 
These were all done by hand calculation. (Note the simi- 
larity with Fig. 5.3 . However, the boundaries at 
a) = 20, 40, 60 were omitted in Fig. 5.3 , so the calcu- 
lations there do not apply to Fig. 6.4 .) It is in- 
teresting to compare the four L n ' s of Fig. 6.4 . The 
larger phase lag of L ^ results in a larger value of 
[dn|L n (0)| -An|L n («)|] , (recall Equation 5.2 and the 
discussion there) . However, since this extra phase lag 
is in the low-frequency range, and the value at to = 10 
is even higher than the rest, this difference means only 
a much larger l L ni(°)l than the rest. One should not 
make much of the small differences at to = 100 rps , be- 
cause the design was performed by hand calculations, 
which do not permit very precise optimum design. 


In any case, the resulting design leads to the 
changes in 


L ( jto) 
1+L ( jto) 


= | T ' (jto) | 


shown in Fig. 6.5a . Most of the responses lie outside 
the permitted boundaries* A prefilter F must be added 
and the requirements on |f| are found, as explained in 
Sec. 6.1 . The results are shown in Fig. 6.5b . One may 
note that over some frequency range (up to about 3.3 rps) 
the boundaries are completely filled, but not for to > 3.3. 
This may seem surprising, because if L(jto) lies pre- 
cisely on its boundary at each frequency, then Aftn|T(jto) | 
should cover precisely its entire permissible range. The 
reason for this discrepancy is that the boundaries have 
been restricted to some extent, due to the constraint on 
the disturbance response damping factor, as explained in 
Sec. 2.3 . 




DEGREES 


Fig. 6.4: Boundaries of L(jw) and some L(jto) designs 
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The step responses resulting from Fig. 6.5(b) are shown 
in Fig. 6.6 . These responses fall within the time response 
bounds corresponding to the complex pole range of Fig. 6.1b. 
However/ there are some time responses whose behavior, as a 
function of time, may not be acceptable. The reference is to 
Curve 8, in which the first overshoot is below the final 
value. The corresponding frequency response marked 1 8 ’ in 

I 

Fig. 6.5(b) has a minimum followed by a peak which does not 
reach the | T (0 ) | reference. We shall refer to this as the 
"wobbling problem". (It is interesting to note that reflec- 
tion of the | T ( j to) | frequency response about the y axis at 
w = 0 , gives a roughly qualitative picture of the system 
step response. This is especially true of the 'wobbling' 
phenomenon.) The wobbling problem will be next considered. 

6 . 3 The Wobbling Problem 

Case 8 in ELg . 6.5b is a good example of a 'wobble' in 
| T ( jco) | . Evidently, at a low frequency, (2 rps) , case 8 

must be at the low end of the template of P ( j 2 ) , see Fig. 

6.3, template marked to = 2 . When this template is posi- 
tioned to find the boundary of acceptable L(j2), Case 8 
lies near the low extreme of j T ( j 2 ) | . On the other hand, 

at to = 3 , if the template of P ( j 3) is calculated, it is 
found that Case 8 lies near the top of the template, towards 
the left side, which means that it will be near the high ex- 
treme of | T ( j 3 ) | . The occurrence of a wobble of this type 

may therefore be predicted when there are some plant condi- 
tions which exhibit this kind of behavior. 

How can the 'wobble' in the frequencu response be eli- 
minated? The simplest but least economical way is to in- 
crease the level of | L | for all w . 


The corrective 
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0 1.0 2.0 


t , sec 


Fig. 6.6: Step responses for T ( j co ) of Fig. 6.5b 



75 . 


effect is as follows — see Fig. 6.7 . At oj = 2 , the 
position of L 8 (j2) is, say, in the neighborhood of Q , 
i.e., in the low end, while at oj = 3 , L 8 (j3) is say, 

at U , in the high end. Increase of |l| at the same 
phase improves matters at both frequencies. It decreases 
| T | at oj = 3 , and increases it at to = 2 , both ef- 
fects helping to straighten out the wobble. This method 
was used in the above design example, increasing j L | by 
2 db . The new |FL/(1+L) | values are shown in Fig. 6.8a, 
and their step responses are sketched in Fig. 6.8b . 

It is certainly possible to eliminate the wobble in a 
more economical manner, i.e., with a smaller increase of 
| L | at high frequencies. The two principal frequencies of 
the wobble, i.e., the minimum and maximum points, which 
were at co = 2,3 for Case 8 in Fig. 6.5b, are considered. 
Fig. 6.9a is used to present the argument, with extreme 
wobbling for Condition 1, between co x and w 2 . Suppose 
L is changed at co x such that |tJ is increased and 
[ T 2 1 is lowered, from points A x , A 2 to points Bj , B 2 , 
respectively. This by itself alleviates the wobbling prob- 
lem in two ways. First, the difference in level between 
T j (ai 1 ) and T 1 (oj 2 ) is reduced by log B : - log A } = 6 x . 

Second, it is possible to improve matters by modifying the 

\ 

prefilter (F) value at w 1 . Let the new | F ( j oj x ) [ be 
larger than the old, by the amount log A 2 - log B 2 = <5 2 • 
The result is to restore | T 2 ( j co x ) j back to the position 
A 2 , and to increase | T x ( j co x ) | precisely by the amount 
S 2 . Thus, the total improvement is 6j + <5 2 . The same 
two effects occur at co 2 if a change of L at oj 2 causes 
a decrease in | T x ( j co 2 ) and an increase in | T 3 ( j a> 2 ) | . 

There remain the following questions to be answered. 
How much of the improvement should be assigned to co 1 , 
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and how much to to 2 ? Whatever the amount assigned to 
either, how is it to be attained, i.e., how should L be 
changed at that frequency? At any frequency one can find 
the net gradient of | T | , as follows. Consider Fig. 6.9b, 

whose points Aj,A 2 correspond to the equivalent points 
in Fig. 6.9a . The gradient of |t| is, of course, normal 
to the contours of constant |t| . The solid arrows are 
drawn in the desired direction in each case, i.e., at A 2 
it is desirable to decrease | T | , while at Aj , it is 
desirable to increase | T A | . Let the length of the arrows 

be assumed proportional to the gradient magnitudes. Then 
the net gradient is the vector sum of the two arrows and 
gives the optimum direction for decreasing J T | with 
least change in | L | , at that specific frequency. Super- 
ficially, it would seem that one could thus find the net 
gradients at each of tOj , to 2 and assign the burden of 
correction to each, in proportion to their net magnitudes. 
But the matter is not that simple, because a cnange in 
| L | is associated to some extent with a change in _/_L 
and the change cannot be local only. For example, suppose 
AL = 2 db / 30° is desired at to < to . This can be 

L 1 

achieved by means of a suitable 'finite line segment' 
properly placed (Fig. 6.9c) . The phase lead exists for 
only a part of the frequency range, while the magnitude 
change persists for all to > to and its benefits are 
therefore available for all to > to . Hence, if the to 

X ^ 

of Fig. 6.9a is not much larger than to , it is cer- 
tainly worthwhile maintaining the level of the correction 
(denoted by G) to to . Also, it is clearly better to 
use Fig. 6.9c , rather than increase the level of |l| by 
the same amount for all to , because Fig. 6.9c permits 
the desired phase lead to be simultaneously obtained, and 
so obtain a correction lined up with the net gradient 
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calculated for the point oi : . It is possible to make addi- 
tional adjustments at a) , to try thereby to achieve a 
correction there, also lined up with the gradient. 

The net 'loss' to the system, in the form of larger 
| L | at high frequencies, is available from Equation (5.2) , 

when | L (0 ) | is maintained at its former value. It is, 
therefore, desirable at higher frequencies where correction 
is not necessary, to reduce |g| back to zero db . If that 
is possible, them there will be an associated phase lag 
which will cancel the phase lead in Fig. 6 . 9 c, and by 
Equation (5.2), there will then be no net loss. This is 
impossible to do completely, if our original design for L 
was indeed an optimum design. But it is at least possible 
to partially do so (see Fig. 6.9d). Let L : be the 

original optimum loop transmission which resulted in the 
wobbling problem. Let the circles represent the result of 
the correction which eliminated the wobbling. For simpli- 
city, it is assumed oj 5 is sufficiently large so that the 
phase of the correction G is back to zero at w 5 , so 
the circles are vertically above the former points on Lj . 

A negative 'finite line-segment' if introduced here, 
would have the effect indicated by the arrows. It could, 
at least theoretically, be shaped optimally to reduce the 
'loss' as much as possible. It is certainly better to do 
this than to allow the correction, i.e., this higher level 
of | L | , to stand as is over this higher frequency range 
where it is not needed. 



CHAPTER 7 


CONCLUSIONS 


7.1 Unfinished Work 


This report has presented a procedure for design of 
two-degree-of-freedom feedback systems with parameter ig- 
norance and specified acceptable time-domain performance 
bounds. One of the attractive features of the technique 
is its transparency. At every point the designer under- 
stands the basic conflicts between the desiderata, espe- 
cially the relation between the loop bandwidth and its 
effect on the system response variations. Another advan- 
tage is in the speed of the design. Once the designer 
understands the method, he can design a problem such as 
the one described in Chapter 6 in about ten working hours, 
if he has computer programs for evaluating time from fre- 
quency response and vice versa. 

There are several important uncompleted parts. One 
is the study of the optimum L(jco) for the general case. 

A second is the wobbling problem. A third is the transla- 
tion of time response bounds to frequency response bounds. 

In addition, it would be very desirable to automate the 
design steps as much as possible. The automation of the 
calculation of the templates of P ( jw) is easy and has 
been done. The determination of the boundaries of L(jco) 
on the Nichols chart is relatively easy to do by hand and 
is probably best left that way. But it would certainly 
be desirable to have a computer program for finding L(jw), 
or at least to polish it up after the initial approxima- 
tion by the methods of Chapter 5. 

7 . 2 Application to Disturbance Attenuation 

The same technique may be applied to the problem of 
disturbance attenuation accompanied by plant parameter ig- 
norance. Consider, the three disturbances of Fig. 7.1, one 
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by one . 



(7.1) 


with 


L = GP ! P 2 


(7.2) 


Let l = j- and then 


C 

D, 


145 , 


(7.3a) 


«.n|^| -tnl^ 


(7.3b) 


The rotation of the Nichols chart by 180° is equivalent 
to the transformation L = 1/i . Hence, the Nichols chart 


may be so used for T. 


Templates of 


p = p 1 p 2 


are found 


in the same manner as before and so are boundaries of i (jto) 
to satisfy given frequency response specifications on 


T , etc. 

U 3 


Consider 


_C_ 
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+ P o G ~ + L o 


where 


P fl is any nominal plant transfer function. It may be 


more convenient to work in the complex plane, rather than in 
the Nichols chart. At any given frequency, the range of. 
P Q /P may be found. Suppose it is as shown in Fig. 7.2 at 


to = to. 


Suppose also that - L 0 (jto 1 ) at point A : is 


under consideration. 


Then 


p o ( 3 W ! ) 
AB 


(7.4) 


where B 


may range anywhere inside or on — ( j to) 


One may 


then find the boundary of the - L (jto ) ' which satisfies the 








frequency response specifications on 
the design is straightforward. 


The balance of 
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Finally, consider 


i o 


2 0 


1 + PjP 2 G 


p p 
^ 1 0 2 0 


p p 
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+ L, 


P 

A io 


The boundaries of Y . = — p— / 

c i 


1 o 


2 0 A 


= Y 


1 2 


are 


1 1 2 

sketched as in’ Fig. 7.3 . The region A of y 1 

the larger region A' of Y 12 • ® into B' and 
A'jB' may overlap. This is all at some specific 


a) . From Equation (7.5) 


‘D. 


^ 2 0 Y ! 

MN 


(7.5) 


obtained and 

maps into 
parts of 
frequency 


(7.6) 


When the region A in is considered, then the appropriate 

range of N , in Fig. 7.3, is A' etc. In this way boundaries 
of acceptable -L Q (jw 1 ) may be obtained and the design con- 
tinued by the methods of this report. Alternatively, one may 
work on the Nichols chart by writing 

P P -i 
12 P L/P, 

m = = L (7.7 

a d 2 1 + P ! P 2 G 1 + l 


so 


l n T t 


= l n 


1+L 


+ l n P 


(7.7b) 


When the template of P is manipulated on the Nichols 
chart, it is necessary to consider also the template of Pj , 
in a manner equivalent to Fig." 7. 3 . 

Finally, if there are specifications on both the system 
response to commands, and its response to disturbances, it is 
possible by the means described to find the resulting 
boundaries of L(joo) for each separately, and then take the 
worst of both, such that both sets of requirements are 
satisfied. 
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